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A circulating current can be induced in the Fermi sea by displacing a scatterer, or more generally 
by integrating a quantum pump into a closed circuit. The induced current may have either the same 
or the opposite sense with respect to the "pushing" direction of the pump. We work out explicit 
expressions for the associated geometric conductance using the Kubo-Dirac monopoles picture, and 
illuminate the connection with the theory of adiabatic passage in multiple path geometry. 



I. INTRODUCTION 



II. OUTLINE 



The question "how much charge is transported due to 
an adiabatic translation of a scatterer" has been raised 
in the context of an open geometry in Refi The scat- 
terer is a potential barrier whose location x = X\ and 
transmission gx = g{X 2 ) at the Fermi energy are deter- 
mined by gate controlled parameters X\ and X 2 . For 
the single mode "wire" of Fig. la, using the Buttiker- 
Thomas-Pretre (BPT) formalism^ 3 -, one obtains the fol- 
lowing result: 



Q = (1 - <? x ) -fe AX i 



(1) 



where fc F = (2m£' F ) 1 / 2 is the Fermi momentum, and 
AX i is the translation distance of the scatterer. Ref. 1 
has referred to this transport mechanism as "snow plow" 
and pointed out that it should be regarded as the pro- 
totype example for quantum pumping: a full pumping 
cycle (Fig. 2a) would consist of translating the scatterer 
to the right, shrinking its "size" , pulling it back to the 
left, and restoring its original size. 

Quantum stirrings^ is the operation of inducing a DC 
circulating current by means of AC periodic driving. This 
is naturally achieved by integrating a quantum pump in a 
closed circuit&££. In particular Refs4£ have considered 
the same adiabatic "snow plow" mechanism as described 
above and obtained for the model system of Fig. lb the 
following result: 
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where g v is the transmission of the ring segment that 
does not include the moving scatterer, as defined by its 
Landauer conductance if it were connected to reservoirs. 

Eq. ((21) is "classical" in the Boltzmann sense because in 
its derivation the interference within the ring is ignored. 
The purpose of the present study is to derive quantum 
results for the stirring in a low dimensional device, where 
quantum mechanics has the most dramatic consequences. 
In particular we would like to illuminate the possibility 
of having a counter-stirring effect: by "pushing" the par- 
ticles (say) anticlockwise, one can induce a circulating 
current in the counter-intuitive (clockwise) direction. 



As a preliminary stage we provide a simple pedagogical 
explanation of the counter-stirring effect by regarding the 
"pushing stage" of the pumping cycle as an adiabatic 
passage in multiple path geometry^. 

For the actual analysis in the general case we use the 
Kubo-Dirac monopoles picture of-£. Within this frame- 
work the pumped charge Q is determined by the flux 
of a B(X) field which is identified as the Berry-Kubo 
curvatur e 11 ' 12 ! 13 . We study both analytically and nu- 
merically how this field looks like. The results are illus- 
trated in Figs. 2-4. Summing the contributions of all the 
occupied levels wc get expressions for the geometric con- 
ductance G. Integrating over a full pumping cycle we get 
results for Q. 

We derive practical estimates for the stirring which is 
induced due to the translation of either small (g x ~ 1) 
or a large (g x <C 1) scatterer, including the possibility 
of having g x ~ <? v . The dependence of Q on the size of 
the scatterer is plotted in Fig. 5, where it is contrasted 
with the classical expectation, and compared with the 
analytical approximations. In the Summary we refer to 
the experimental measurement aspect. 



III. THE COUNTER-STIRRING EFFECT 

The essence of the counter-stirring effect can be un- 
derstood without the Kubo-Dirac monopoles picture by 
adopting the "splitting ratio" concept of Ref.- . Referring 
to Fig. lb the translation of the scatterer to the right is ef- 
fectively like lowering the potential floor in the left bond, 
and raising the potential floor in the right bond. This 
induces an adiabatic passage of a particle from the right 
to the left. The particle has two possible ways to make 
the passage: either via the "V" barrier (coupling W^ 2 ) or 
via the "X" barrier (coupling Wf 2 ). The splitting ratio 
determines the fraction of the current that goes via the 
"V" barrier: 
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where the last equality is based on the analysis in Ref£. 
If W\ 2 were the classical rate of the transition, we would 
have < A < 1 and the current would flow in accordance 
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with our classical intuition. But W^ 2 and Wf 2 ar e real 
amplitudes that might have opposite signs if an odd level 
crosses an even level. Consequently if \Wf 2 \ > I W12 1 we 
get A < which implies that a circulating current is in- 
duced in the counter-intuitive (clockwise) direction. This 
does not come in any contradiction with the observation 
that the net transport (summing over both barriers) is 
still from right to left. 



IV. THE MODEL HAMILTONIAN 

Our model is a ID coherent ring with a fixed scat- 
terer and a controlled scatterer. The fixed scatterer is 
some potential barrier V(x), and the controlled scatterer 
is modeled as a delta function whose position and trans- 
mission are determined by the control parameters X\ and 
X 2 respectively. The one particle Hamiltonian is: 



H 



±-f + V(x) 
2m 



■ X 2 (t)5(x - X x (t)) 



(4) 



with periodic boundary conditions over x G [— L/2,L/2] 
so as to have a ring geometry. Below we fur- 
ther assume that both bonds are of similar length 
(L— X\ ~ X% ~ L/2). The current is measured through 
a section x=xq= + at the fixed barrier, and accordingly: 



1 = — (pS(x - x ) + 5(x - x )p) 
zm 



(5) 



We also define generalized forces which are associated 
with the control parameters: 
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For practical use it is more convenient to describe the 
fixed scatterer by its scattering matrix, which can be 
written as: 



S v = 
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(8) 



We study the case when the model parameters are such 
that the transmission of the fixed scatterer is small 
(<7v 1) 5 the Fermi momentum is large (k F L 3> 1), and 
the controlled scatterer is translated a distance AXi that 
equals several Fermi wavelengths. 



where G 31 is called the geometric conductance. For two 
parameters driving one can write: 
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where G = (G 31 ,G 32 ), and X = (X U X 2 ) and ds = 
(dX 2 , — dXi). For a particle that evolves adiabatically in 
the level n we have G 31 = B 2 and G 32 = —B\ where: 
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In fact (Bi,B 2 ) are elements of the Kubo-Berry 
curvatur e 11 ' 12 ! 13 which one can regard as a fictitious 
magnetic field B = (Bi,B 2 ,B 3 ) in an embedding space 
X = (Xi, X 2 , X 3 ). From the requirement of having well 
defined Berry phase it follows that the sources of B(X), 
which are located at points of degeneracy are quantized, 
so called "Dirac monopoles" . 



VI. THE X SPACE 

Due to the gauge symmetry $ $ + 2irh/e the Dirac 
monopoles are arranged as vertical chains (see Fig. 2(f)) 
[7i=e=l]. Due to the time reversal invariance of our 
7i{Xi : X2) it follows that a Dirac chain is either a du- 
plication of in-plane monopole at X3 — or off-plane 
monopole at X3 = tt. Let us find an explicit formula 
for the (Xi,X 2 ) locations of these vertical chains. The 
equation for the adiabatic energies E n (X) is of the form 
cos(fcL+7) = v /gcos($) where g is the total transmission 
of the ring and 7 is the total phase shifts of the scatterers 
(the fixed scatterers plus the moving scatterer). If X 2 is 
such that g x (E; X 2 ) — g v (E) we can always find X\ such 
that the total transmission would be g — 1, which is the 
necessary condition for having a degeneracy. Together 
with the equation k E L + "/(E) = m with r=integer this 
defines a set of energies E r — (k r ) 2 /2m and associated 
values X 2 for which the n=r level has degeneracy with 
the n=r+l level provided X\ is adjusted. To be more 
precise r=even are in-plane ($=0) degeneracies, while 
r=odd are off-plane (<E>=7r) degeneracies. The X\ loca- 
tions of these degeneracies are half Dc-Broglie wavelength 
apart (see Fig. 3): 
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If we were changing the flux X3 = $ through the 
ring, the induced current would be given by Ohm law 
I = (I) = — G 33 Xs, where G 33 is the Ohmic conductance 
and — X3 is the electro-motive-force. Similarly for a vari- 
ation of the parameter X\, the current is I = — G 31 Xi, 



where the [1/2] shift applies to in-plane degeneracies. 
The arrangement of the degeneracies in X space is il- 
lustrated in Fig. 2. For each (X"[ : X 2 ) we have a vertical 
Dirac chain whose monopoles are formally like sources 
for the B field. 
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VII. FERMI OCCUPATION 



X. MATRIX ELEMENTS 



If we have many body system of N = ^ n f n particles, 
then B = f n B^ n \ At finite temperature each occu- 
pied level (except n=l) contributes two sets of X£) 
chains, namely r—n and r=n— 1, which are associated 
with the E n = E n ±\ degeneracies. By inspection of 
Eq. tfTU)) . taking into account that T nm is antisymmet- 
ric, one observes that the net contribution of the rth set 
of Dirac chains is f r — f r +i- In particular for zero tem- 
perature Fermi occupation, the net contributions comes 
from only one set of chains which is associated with the 
degeneracies of the last occupied level with the first non- 
occupied level (Fig. 2bcd). 

VIII. CLASSICAL LIMIT 

At finite temperatures we can define the smeared prob- 
ability distribution of the Dirac monopoles with respect 
to X 2 as follows: 



f(X 2 ) = J2lfr-fr+l}S(X 2 -Xr) 



(12) 



Disregarding fluctuations Eq. (J9j> implies a monotonic de- 
pendence of Q on X2 in qualitative agreement with 
Eq. ([2]) . If the expression in the square brackets of Eq. 

were equal J Q 2 f(X')dX' , it would imply a quantitative 
agreement as well. In order to have this quantitative 
agreement we have to further assume that the distribu- 
tion f(X 2 ) is determined by some chaotic dynamics in 
the scattering region which would imply erratic depen- 
dence of the S matrix on the energy E. See^ for further 
discussion. 



IX. QUANTUM LIMIT 

Our interest below is in the opposite limit of zero tem- 
perature where f(X 2 ) becomes a step function. Ob- 
viously in this limit a step like behavior of Q versus 
X 2 would be a crude approximation. By inspection of 
Ea.tflO)) it follows that the result for G = G 31 = B 2 is 
very well approximated by 
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where n is the last occupied level. This observation, as 
well as the associated analytical results which are based 
on it, have been verified against the exact numerical re- 
sults of Figs. 3-4. Below we derive explicit expressions 
for G vs Xi for both small and large values of X 2 . Our 
results for Q are plotted in Fig. 5. Note that the depen- 
dence on X\ has ir/k F periodicity due to the X space 
arrangement of the monopoles, and accordingly the inte- 
gration gives Q tx ek F AXi/ir, with a prefactor that we 
would like to estimate. 



The matrix elements of the current operator X and of 
the generalized force T x are 



X'nr, 
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(14) 



(15) 



where dtp = l/2(dipL + dipn) is the average derivative on 
the left and right sides of the delta barrier, and the second 
expression for T^ m was obtained by using the matching 
conditions across the delta barrier. The wave function is 
written as ip{x) = Csin(<^ + kx). We found that a very 
good approximation for T nm with m=n+l is 



(16) 



where the + (— ) sign is for n =even (odd), and v F = 
fc F /m is the velocity in the energy range of interest. For 
the calculation of Tl nn and E m — E n we have to distin- 
guish between the two cases of small/large scatter. This 
means the small/large X 2 regimes where g x ~ 1 or <C 1. 



XI. TRANSLATING A SMALL SCATTERER 

If the controlled scatterer is small, we treat it as a 
perturbation. For the energy splitting we get 

E m -E n « ^v F T \x 2 cos(2k F X 1 ) (17) 

where for notational convenience we take X[ as the new 
origin. After some further algebra we get 

2k 7T 
T x mn = ±X 2 ^- cos (2fc F Ax) +X 2 jz (18) 

where the ± sign is as in Eq. (|16p . The conductance can 
be written as 



G = —k F \^ Q v cos [y 2k F X\) 



(19) 



where the coefficients of the leading non-negligible terms 
[the small parameter being (1— <? x )/<?x] are 



Go = ±2yfgZ 
Gx 
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(20) 
(21) 



Upon integration we get Q = —&Qq per half Fermi wave- 
length displacement of the scatterer. 
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XII. TRANSLATING A LARGE SCATTERER 

If the controlled scatterer is large, most of the charge 
transfer is induced during the avoided crossings (sharp 
peaks in Fig. 3 lower panel). Consequently we use the two 
level approximation scheme ofS. with m = n+1 leading to 
the results: 



E m — E„ 
T 1 



\R\ 



2 2-^2 

±1™* W\ 



(22) 
(23) 



where the ± sign is as in Eq. lfTB)) and the dimcnsionlcss 
distance in X space from the degeneracy point is: 



R = \2k F {X x - XI), 

Accordingly the conductance is 

fc F R 2 



*9\ 



\{X 2 



G = 



(24) 



(25) 



Integrating over X\ we get Q — e\(X2) per half Fermi 
wavelength displacement of the scatterer, as expected 
from the "splitting ratio" argument. 



consequences would be observed in low dimensional de- 
vices. For this reason we have analyzed in this paper the 
prototype problem of "pushing" a current by translat- 
ing a scatterer in a single mode wire. We have obtained 
explicit results for the B field, which determines the geo- 
metric conductance G, and consequently the Q of a closed 
pumping cycle. We also illuminated the counter-stirring 
effect using the splitting ratio concept of adiabatic pas- 
sage in multiple path geometry. 

A few words are in order regarding the measurement 
procedure and the experimental relevance. It should be 
clear that to measure current in a closed circuit requires 
special technique s 14 ' 15 ! 16 . These techniques are typically 
used in order to measure persistent currents, which are 
zero order (conservative) effect, while in the present pa- 
per we were discussing driven currents, which are a first- 
order (geometric) effect. It is of course also possible to 
measure the dissipative conductance (as in 14 ). During 
the measurement the coupling to the system should be 
small. These are so called weak measurement conditions. 
More ambitious would be to measure the counting statis- 
tics, i.e. also the second moment of Q as discussed 
which is completely analogous to the discussion of noise 
measurements in open system s 18 ! 19 . Finally it should be 
pointed out that the formalism above, and hence the re- 
sults, might apply to experiments with superconducting 
circuits (see£). 



XIII. SUMMARY 



The integration of a two-terminal quantum pump in a 
closed circuit is not a straightforward procedure. Due to 
interference the pumped charge Q would not be the same 
as in the Landauer/BPT setup, and even the sense of the 
induced current might be reversed. The most dramatic 
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FIG. 1: In the case of an open geometry the pumping device is 
connected between two unbiased reservoirs (panel (a)), while 
in the present study it is integrated into a ring (panel (b)). 
The induced current is measured through a section indicated 
by a dashed line. See the text for further details. 
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FIG. 2: (a) Three representative pumping cycles. Several 
sets of Dirac monopoles may have non-zero weight depending 
on the occupation. Panel (b) is for single level occupation 
where two sets have non zero weight, while either (c) or (d) 
are for zero temperature Fermi occupation. Filled (hollow) 
circles indicate in(off )-plane monopoles. Panels (e-f) give a 
detailed illustration of the associated B field, as implied by 
the numerical findings of Figs. 3-4. The X\ tick marks in (e) 
are half De-Broglie spaced, while the horizontal blue lines are 
paths for which numerical results are presented in Figs. 3-4. 
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FIG. 3: The conductance G=G 31 of Eq[T0]is numerically cal- 
culated for a ring of length L=151.43 with V(x)—U5(x) where 
t/~10. We consider single level occupation n=138. At this 
energy g v —0.06. The upper (lower) panel is for translation 
of a very small (large) scatterer with g x =0.98 (gx~8 ■ 10~ 8 ) 
corresponding to the lower (upper) horizontal blue paths that 
are indicated in Fig. 2e (same X\ axis). 
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FIG. 4: Additional plots of the conductance G as calculated 
in the previous figure. The upper (lower) panel is for trans- 
lation of a scatterer with g x =0.20 (g x =0.03) corresponding 
to the horizontal blue paths in Fig. 2e that go below (above) 
the <7x=gv axis. Note that the large peaks are positive (neg- 
ative) while the small positive peaks switch sign only when 
the scatterer is lowered further. This indicates that the field 
lines bend in the X% direction as illustrated in Fig.2f. 
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FIG. 5: Numerical calculation of Q for zero temperature 
Fermi occupation: all the levels are occupied up to n=137 
(upper panel), and up to n=138 (lower regular and zoomed 
panels). The model parameters are the same as in Fig. 3. 
The integration was carried out along segments similar to the 
paths shown in Fig. 2e where g x varies between (? X =8.30T0 -5 
and (/ x =0.45. For sake of comparison we display both the 
analytical classical (Eq[2]) and quantum results (Eqs l20& i l25|) . 
The value of X2 for which <?x=<?v is indicated by a vertical 
line. 



